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It is shown that the distribution derived from the principle of maximum Tsallis entropy is a 
superposable Levy-stable distribution. Concomitantly, the leading order correction to the limit 
distribution is also deduced. This demonstration fills an important gap in the derivation of the 
Levy-stable distribution from the nonextensive statistical framework. 
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Recently, Prato and Tsallis [1] gave an argument to demonstrate the appearance of the Levy- 
like distribution function from the nonextensive statistical framework applied to the problems of 
random walk and diffusion phenomena. The distribution they obtained has the same asymptotic 
behavior as that of the Levy distribution for large values of the random variable. Then, 
employing the Levy-Gnedenko generalized central limit theorem [2,3], they claim that their 
distribution approaches the exact Levy distribution by convolution of the distribution many 
times. However, an important question remains regarding the convergence of the convolution. 
The purpose of this paper is to examine this question and to explore the rate of convergence. 

In Ref. [1], the authors consider the maximization of the Tsallis entropy [4] 

s M = ^- q {] f (i) 



under the constraints on the normalization condition of the distribution 



| dx p(x) = 1 , (2) 



and on the generalized second moment [5] 



dxx 2 P{x) = o\ (3) 



where a is a positive constant giving a length scale and P(x) is the escort distribution associated 
with p(x) [6] 



PM= . [P(X)Y ■ (4) 



\ dx'[p{x)] q 
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q in eqs. (1) and (4) is a parameter, which is assumed here to be 



<q<3. 



(5) 



The maximum-entropy distribution is found to be 



p(x) 



1 4-1 



3--? If, . g-1 x 2 ^ 
3-q O 2 



ylq-2j 



1 + 



1/(9-1) 



(6) 



The range in eq. (5) covers the same range as that of the Levy distribution. [See eq. (10) below.] 
It should be noted that the ordinary second moment with respect to this distribution is divergent 
but the generalized second moment in eq. (3) is finite. For large values of x, the above 
distribution asymptotically behaves as 



p{x) ~ x~ 2l{q - l) . 



(7) 



Comparing this with the symmetric Levy distribution 



L r (x) = — | dk exp(-ikx) exp(-a| k\ 7 } (0</<2) 

2 7L 



(8) 



with a positive constant a , whose asymptotic behavior is 



L(x)~x- l -\ 



(9) 



we find that q is related to the Levy index J as follows: 
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q = 



3+ y 
1+7 



(10) 



The generalized central limit theorem is stated as follows. Given a distribution whose 

ordinary second moment is divergent, its iV-fold convolution approaches one of a stable class of 
distributions r (x): < y < for large ./V provided that the convolution converges. This 

property has not been ascertained in Ref. [1]. In what follows, we discuss this problem using the 

method of characteristic functions. 

In order to apply the generalized central limit theorem, it is necessary to consider the sum of 
N scaled independent random variables {*,} i=1 2 L N [2,3]: 



x 



(11) 



The scaling factor B N here has to be chosen in such a way that the limit distribution is 
independent of the number N of convolutions. Suppose each x t obey the Levy distribution in eq. 
(8). Then, the distribution of x is given by 



L^(x) = B N (L y *L *L y ){B n x). 



(12) 



The characteristic function of L ( Y } (x) is calculated to be 



X[ N) (k)= ] dxcxp{ikx)Lf(x) 









N 


Xl 


\ B Nj 













(13) 



where % L {k) is the characteristic function of L y (x): 
X L (k) = exp(-a\k\ 7 ). 



(14) 
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Therefore, we have 



B N = N 



(15) 



where the constant of proportionality is set equal to unity for the sake of simplicity. 

To examine the convergence properties of the distribution in eq. (6), first we calculate its 
characteristic function. After some algebra, we find 



X(k) = jjjc exp(ikx) p(x) 



(16) 



where T(z) and K v (z) are the Euler gamma function and the modified Bessel function [7], 
respectively, and 



3-q y 



2q-2 2 



(17) 



A = a/2vct. 



(18) 



Clearly, the range of v is 



0< v<l. 



(19) 



Now, let us discuss the limit /V — > oo of 



/(*; n) - 



x 



(20) 



in order to determine the convergence properties of the distribution in eq. (6). For this purpose, 
we take the logarithm of f(k; N): 



lnf(k; N) = N In 



, l-y/2 



r k\k\ 

N llr , 
\i\ j 



r 









K 



r/2 



(21) 



Using the formula K v (z) = — — v -^- ( — 4-^ [7], where I v (z) is essentially the Bessel function of 



2 sin(v7r) 

imaginary argument, and evaluating the series expansion I v (z) = ^ 



{zl 2) 



2k+v 



t k\r(k+ v+i) 



[7], we find 



In f(k; N) = N In 



1- 



1 



1- 



7 



/V rfi + |h 



+ 



i 



r x\k\^ 2 



2-yN 



21 y 



V 2 j 



+L 



7 



V 



fi + i 

I 2 



v 2 , 



+ 



l-yN 2 ' 7 - 1 



v 2 , 



1 

2N 



1-1 



v 



-| 2 



f 7 
1 + ^ 

I 2 



AJ/k 
• 2 , 



+L 



(22) 



Therefore, in the limit iV — > °°, we obtain the exact Levy distribution with 



1-| 



fl + l 

V 2 



f-Y 

v2j ' 



(23) 
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Simultaneously, we also find that dominant correction to the Levy distribution is of O^N l ) if 
< 7 < 1 whereas it is of 0(N l - 2 ' Y ) if 1 < J < 2. 

In conclusion, we have shown that the distribution derived from the principle of maximum 
Tsallis entropy is in fact a superposable Levy-stable distribution. We have also determined the 
leading order correction to this limit distribution. We have thus filled an important gap in the 
derivation of the Levy-stable distribution from the nonextensive statistical approach. 
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